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Abstract 

In extensions of the Standard Model with SU(2) singlet scalar helds, there can be regions of 
parameter space for which the electroweak phase transition is hrst order already at the mean-field 
level of analysis. We show that in this case the phase interface (bubble wall) can become ultra- 
relativistic, with the relativistic gamma factor 7 = (1 — ^aii)" 1 ^ 2 g row i n g linearly with the wall's 
propagation distance. We provide a simple criterion for determining whether the bubble wall "runs 
away" in this way or if 7 approaches a terminal value. 



I. INTRODUCTION 



Electroweak symmetry breaking is a central feature of the Standard Model, arising be- 
cause the Lagrangian mass squared for the Higgs field is negative. However, interaction with 
a thermal medium rather generically raises the effective (thermal) mass squared of scalar 
particles, leading to high temperature symmetry restoration [l|, 0] . Such symmetry restora- 
tion is expected in the Standard Model at temperatures of order T ~ m H ~ lOOGeV, which 
presumably occurred very early in the hot Big Bang. 

Such symmetry breaking in a gauge-singlet scalar field would automatically imply a ther- 
mal phase transition between symmetry-restored (high temperature) and symmetry broken 
(low temperature) phases. 1 But the Higgs field (or fields) transforms nontrivially under a 
nonabelian gauge symmetry. The notions of vacuum expectation value and symmetry break- 
ing for such a scalar are only defined perturbatively, and there is no guarantee that there 
will be a phase transition associated with electroweak symmetry breaking [i| . In particular, 
within the Minimal Standard Model at allowed values for the Higgs mass, no electroweak 
phase transition exists [5|. But an electroweak phase transition would occur in numerous 
phenomenologically interesting extensions of the Standard Model. Such a transition is in- 
teresting cosmologically for a number of reasons: 

• A first order electroweak phase transition might lead to electroweak baryogenesis jf], 0] , 
giving rise to the baryon number observed in the universe today. 

• A first order phase transition might produce primordial magnetic fields jil, 0] which 
might seed the formation of today's galactic and extragalactic magnetic fields. 

• A first order electrowe ak p hase transition might give rise to observable gravitational 
wave backgrounds [lo|, [lj [l2| in interesting frequency ranges for future experiments 
(around 10~ 3 Hertz). 

There are two key attributes of the transition which are important in these scenarios. One 
is the strength of the electroweak phase transition (for instance, the available free energy 
and latent heat as a fraction of the energy density of the electroweak plasma). This is 
determined by model parameters (field content, masses and coupling constants) which are 
not known. However we have developed the theoretical tools for determining the strength 



of the transition, given the model parameters [13(. The other attribute is the bubble wall 
velocity, meaning the propagation speed of the interface between high temperature and low 
temperature phases, while the transition is taking place. This velocity is relevant to all 
physical consequences of the phase transition. In particular: 

• Electroweak baryogenesis scenarios seem to be quite sensitive to the bubble wall ve- 



locity [14], ll5|, UM, \Y7\, |18|]. A slow bubble wall allows nonequilibrium conditions to 



exist "in front" of the bubble wall, in the symmetric electroweak phase where baryon 



number violating processes are fast [14| . A thin, very fast bubble wall would provide a 



1 It is possible in some models to have thermal effects induce symmetry breaking, rather than restoration 
[jj. But this requires rather special coupling relations and is not important to the rest of our discussion. 
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nonequilibrium environment only in the broken phase, probably making baryogenesis 
very inefficient. 



Gravitational waves are of maximum strength if the bubble walls are highly relativis- 
tic. Slow subsonic bubble walls give rise to rather extended and smooth hydrodynamic 
waves, which should be poor sources of gravitational waves (see jl9[ for the hydro- 
dynamics of transitions with finite bubble wall velocities). On the other hand, if the 
bubble walls propagate ultra-relativistically and a nontrivial fraction of the free energy 
density of the plasma goes into their acceleration, it could give rise to very interesting 
gravitational wave signals. 2 



In a first-order cosmological phase transition in vacuum [22j, the energy density of the false 
vacuum is completely converted into kinetic energy of the bubble walls. As a result, the wall 
velocity keeps increasing and approaches the speed of light (7 = 1/ y/1 — v 2 grows linearly 
with the distance a bubble wall propagates). Either this behavior, or highly relativistic walls 
but with a terminal 7 factor, have often been assumed in the gravitational wave literature 



11]. On the contrary, the electroweak baryogenesis literature has generally assumed that 



interactions with the plasma impede the bubble wall and keep its velocity subsonic, y < 



1 A/3, as the rather limited literature computing the wall velocity suggests v < 0.2 23|, [24, 



25|, [26|, 123, |2& [29 



The goal of this paper is to determine under what situations the electroweak bubble wall 
can "run away" and propagate with macroscopic gamma factor 7 = (1 — v 2 )^ 1 ^ 2 ^> 1, and 
when it more modestly propagates with jv < 1. At the thermodynamic level, we show 
in Appendix [S] that runaway walls are logically possible under quite weak conditions; but 
whether they actually occur depends on how much entropy is produced by the passing wall, 
which is a problem in the microscopic wall dynamics. However, the relevant physics is 
actually quite simple, and the right tools already exist in the literature. This paper could 
have been written more than 10 years ago. However we find that the results neither appear 
in the literature nor are well known by either community. 

Generally, at the mean field level models without electroweak SU(2) singlet scalar fields 
predict a second order electroweak phase transition. The transition in such a theory may be 
a so-called "fluctuation induced" first order transition (as for the Standard Model with an 
unphysically light Higgs field). We will show that the bubble wall never runs away in such 
transitions. On the other hand, theories with SU(2) singlet scalars can in some cases provide 
a first-order electroweak (SU(2) breaking) phase transition already at the mean-field level of 
analysis. This is true for instance in the nMSSM j3o| (though not for the MSSM). Indeed, 
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Steinhardt has argued [20( that bubble walls cannot approach the speed of light if there is a plasma 
present. However this result was based on assuming that the Chapman- Jouguet condition is fulfilled, that 
is, that the plasma behind the wall moves at the speed of sound with respect to the wall. This assumption 
makes sense for burning detonations, where the (chemical) free energy is released behind rather than at the 
shock, in a region thick enough to apply hydrodynamics; then the fluid must be subsonic in order to push 



the shock. But, as carefully demonstrated by Laine 2l| , this assumption is not valid for cosmological phase 



transitions, where the free energy is liberated internally in the wall; as a result more general detonation 
solutions are possible. We complete Laine's treatment in the Appendix, where we show that runaway 
walls are thermodynamically consistent under quite weak conditions. 
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as we discuss below, it is true already for the simplest extension of the Standard Model, a 
theory with an additional real singlet scalar. We will show the such "mean-field first order" 
electroweak transitions generically, but not invariably, have "runaway" bubble walls, with 7 
growing linearly with wall propagation distance. We present a simple criterion to determine 
whether or not such "runaway" will occur. 

Specifically, the phase structure is determined by an effective potential V(h) which is the 
sum of a vacuum part V vac (h) and a thermal contribution V T (h). The transition occurs when 
V(h) = V vac + V T has two local minima, a shallower "symmetric" s-minimum and a deeper 
"broken" /i-minimum; the exact criterion for how much deeper the h minimum need be is 
determined by a bubble nucleation calculation (3l|. If replacing V T (h) with its 2-nd order 
Taylor approximation 3 (h — h sjm ) 2 V" (h sym )/2 removes the "broken" minimum or raises it 
above the symmetric one, the bubble wall cannot run away. If, after this replacement, the 
"broken" minimum remains deeper than the symmetric one, then the bubble wall will run 
away. The rest of the paper derives and clarifies this result. 

The paper is organized as follows. In Section [III we review briefly how the 1-loop thermal 
effective potential arises. Section II III gives a simple derivation of the friction on the elec- 
troweak bubble wall when it is propagating ultra-relativistically with respect to the plasma. 
Section [TV] illustrates our result using a specific model (a toy version of the nMSSM), showing 
how to determine whether or not runaway occurs as a function of model parameters. Section 
[V] presents some conclusions. There is also an appendix, which discusses the thermodynamic 
criteria under which runaway walls are consistent. 



II. THERMAL EFFECTIVE POTENTIAL: A REVIEW 

Our arguments will be based on the way that the thermal effective potential arises from 
the particles making up the thermal bath. Therefore we start by reviewing this physics. 
Consider for simplicity a theory of two real scalar fields h, s with tree-level effective potential 

V„(h, s) = ^h 2 + ^s 2 + ^h* + ^ + \V . (2.1) 

At zero temperature we determine the vacuum configuration by seeking the global minimum 
of the potential with respect to h, s. If, for instance, m 2 H is negative but m 2 and A are positive, 
the global minimum will be at ho = \J — m 2 H /\ H and Sq = 0. In this toy model such a nonzero 
value spontaneously breaks the Z 2 x Z 2 symmetry group down to Z 2 . 

At finite temperature we must take into account, besides the energy density arising from 
Vva.c(h, s), also the energy density of thermal excitations of these (and any other) fields. 
Fortunately one may define a thermal contribution to the potential V T (h, s) describing the 
free energy density of the plasma. The total potential V = V vac + V T describes the total free 
energy as a function of h, s, and its minimum describes the most favored classical background. 
One efficient way to compute the thermal contributions to V T (h,s) is differentially (the 
tadpole method [3(); we compute the derivatives dV T (h,s)/dh and dV T (h, s)/ds, which can 



3 When multiple scalar fields change VEVs at the transition, this expression generalizes to the quadratic 
form I Yljj SiSjCpV/dsidsj , where Sj runs over all scalars which change VEVs. 
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be integrated to recover the full potential. In particular, to compute the difference between 
V T (h, s) and V T (h+Sh, s), calculate the change in all particle energies in changing the value 
oih: 

m f d3p f(F \ dE ^ a Sh - V" f d3p f^ E p^) dm 2 

5Vt -\] T^ f ^ E ^— 5h -\<j w-w^T-dh Sh > (2 - 2) 

where the sum ^ a is over all fields (for us, just h, s), E P: h, a = \/v 2 + m l(h, s) and m 2 (h, s) 
is the h, s dependent mass squared for particle type a. For a scalar field m 2 is the second 
derivative of V(h, s) with respect to the scalar field value; for instance, in our toy model 
m? s (h, s) = m\ + 3A s s 2 + Xh 2 /2. 

In addition to this energy cost, there is a free energy cost associated with changing the 
particle occupancies f B {E p ^ ta ) from the values relevant at field value h to those relevant 
at field value h+Sh. But this free energy cost is quadratic in Sh, since the equilibrium 
occupancies are a minimum of the free energy. 4 Therefore, if we treat 5h as infinitesimal we 
may drop this term, and reconstruct the thermal contribution to the effective potential by 
integrating Eq. (12. 2p - provided we account for the way that E P)h>a and f B (E p ^ a ) vary with 
h as we integrate. 

Often the particle masses m a (h, s) are small compared to 7rT. Then we can approximate 
f B (E)/2E in the second expression in Eq. (12.21) with f B (p)/2p. In this case we can perform 
the integral: 

dVr = y^ dmZ f d 3 p f B (E pAa ) _ x ^ dm| f d 3 p f B (p) = x ^ dm^T^ 
dh ^ dh J (2tt) 3 2E pAa ~ ^ dh J (2tt) 3 2p ^ dh 24 1 ' ' 

which is trivial to integrate: 



a 



with Vq(T) an undetermined (and for us irrelevant) integration constant. A slightly less 
drastic approximation, useful if m 2 (hi) — m 2 (h,2) T 2 but valid even if m 2 itself is relatively 
large, is to neglect the h dependence of the particle masses in the expression for E, replacing 
E Pi h, a with E Pi h lja - 



v T {h 2 )-v T {h x )^Y,{<^)-<^)) I 



d 3 p f B (E pAua ) 



(2vr)3 2E pMr 



This approximation is the same as determining V T (h2) by performing a 2'nd order Taylor 
series expansion of V T about hi - provided that the potential is at most quartic, and particle 
squared masses quadratic, in the condensate values of scalar fields (true at tree level in 
renormalizable theories). In what follows we will call this the mean field approximation to 
the effective potential. In our toy model, if only the h field takes a nonzero background 
value, the effective mass for the s field is d 2 V vac /ds 2 = m 2 + Xh 2 /2. In this case the s 



4 This is an example of the usual perturbation theory result that changes to the state first give rise to 
shifts in the (free) energy at second order in a perturbation. 
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FIG. 1: Example of the potential V — V vac + V T in our toy model, for parameters which give a first order 
transition. The curves are V(h) at a series of temperatures from high (top) to low (bottom); the dotted 
curves are the potential when the h — phase becomes spinodally unstable (To) and where the h ^ phase 
becomes unstable (Ti). The solid curve is the potential at the equilibrium temperature (T c ) where both 
phases have the same minimum value. 

field part of the thermal contribution to the potential, making the stronger approximation 
Eq. (E3D, would be V T (h) = T 2 m|/24 + T 2 Xh 2 /A8 (plus a constant). The first term here 
is also a constant; the second looks like a correction (positive if A > 0) to the h field mass 
squared. At temperatures where T 2 A/24 > —m 2 H) there is a (meta)stable minimum at h — 
and symmetry will be restored. 

Often this mean-field approximation is insufficient. For instance, in this toy example it 
predicts that the phase transition should be of second order. Continuing to treat m a <C nT, 
but being more careful to account for the difference between f B (E)/E and f B (p)/p, one finds 
a series expansion in m a /7iT: 



which can predict a first-order phase transition, as illustrated in Figure [U Since this transi- 
tion is only first order because of changes in f B {E) / E (changes in the size of the fluctuations), 
it is called a fluctuation induced phase transition. 

There are two potentially significant corrections to this simplistic description. One is 
that it has treated thermal particles as free, neglecting their mutual interactions. These give 
rise to higher order ( "two loop" ) corrections to the thermal potential computed here. These 
corrections can be important in establishing the order and strength of a phase transition 
when it is second or weakly first order. Second, there are subtleties in defining an effective 
potential if a scalar which may take on a vacuum value is a non-singlet under gauge inter- 
actions. Fortunately, for the electroweak transition to be phenomenologically interesting, 
it must be relatively strong. In this case neither of these subtleties is very important and 




(2.6) 
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FIG. 2: Cartoon of a bubble wall: the scalar field ft, as a function of position h(z). Particles "hit" the wall 
from either side, inducing net forces. 



the one-loop treatment described here is generally adequate (unless the vacuum theory is 
strongly coupled, in which case we have little to say). 



III. FRICTION ON A RELATIVISTIC BUBBLE WALL 

When an effective potential has two local minima, there are two possible phases and 
there can be a phase interface between them, where the field varies from one to the other 
value, illustrated in Figure EJ Plasma particles approach the interface from either side. If 
the interface is not too sharp, particles follow classical trajectories in the wall background 
(wave packets are well described by the WKB approximation). In the interface rest frame, 
a particle's energy and transverse momentum are conserved, so the classical trajectory is 
given by 

dr p dE 2 2p z dp z dm 2 dm 2 dh dz 



dt E ' dt dt dt dh dz dt 

But dz/dt = Pz/E, so 



dp z _ 1 dm 2 dh 
dt ~ 2E dh dz' ^ ' ' 

Note that dp z /dt represents a net force on a plasma particle. The total force on all plasma 
particles, per unit surface area of the wall, is 

with f a (p, z) the occupancy of species a with momentum p at distance z from the wall, a may 
be bosonic or fermionic and f a is in general out of equilibrium. This is the force on particles 
due to the scalar background; there is an equal and opposite force on the background, which 
acts as a pressure difference restraining the motion of the interface. The wall also feels a 
force due to the vacuum potential; F va _ c /A = —AV va _ c , pushing the wall forward. Therefore 
the total force pushing the phase interface forward is (V vac (hi) — V vac (h2)) — F/A. 
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In general the occupancies f a (p, z) are all out of equilibrium and it is a nontrivial problem 
to determine them, simultaneously with the interface velocity, shape, and tension. This is 



what previous literature has attempted to do 23|, |24], |25|, |26j, |27|, |28|, |29| . However, two limits 
are very simple. 



The first simple limit is a bubble wall at rest in the plasma frame [26j. In this case, 
energy conservation, together with Liouville's theorem, imply that f a remains always in 
equilibrium. Therefore Eq. (I3.3P becomes 

with f a the equilibrium distribution function. This expression coincides with the expression 
we found earlier for the thermal potential difference between phases V T {h,2) — V T (hi), so 
the total force on the wall is V(h{) — V(h?). The critical temperature T c is defined as the 
temperature where the potential is equal in the two minima. Therefore, at T c the forces are 
in balance and a wall at rest will remain at rest, as expected. 

The other simple limit is the limit we need in this paper; a wall with velocity v which is 
ultra-relativistic in the plasma frame, i.e., with gamma factor 7 > 1. In practice we will be 
interested in really macroscopic values 7 ~ 10 9 , so in the following we make no qualms about 
strictly expanding to lowest nontrivial order in 7. We begin by analyzing what happens in 
the rest-frame of the bubble wall. In this frame the plasma is ultra-relativistic with velocity 
v, which we take to point along the — z direction (so the plasma-frame wall has positive v z 
velocity). The mean wall- frame particle energy is E ~ 7T, which will always be enough 5 to 
allow the particle to pass over the wall, 7T ^> m a (h 2 ). Furthermore, the mean momentum 
p z ~ 7T is so large that the WKB approximation is excellent and the reflection coefficients 
are exponentially suppressed. Interactions or scatterings between plasma particles occur 
at a time-dilated rate t ~ I/7T (times inverse powers of couplings) and can be neglected. 
Therefore the occupancies evolve undisturbed. Also, the number of particles approaching 
the wall from the positive z side is exponentially suppressed by exp(— jm a /T). That means 
that the particles approaching the wall have received no signal that the wall is approaching, 
and are in equilibrium. 

Rather than integrate the pressure per length over the wall, we can directly compute the 
net pressure on the wall by summing up the momentum change of each particle which passes 



over the wall 26 



j = j (^ /a(p ' in) ( Pz ' in ~ Pz ' out ) ■ (3 - 5) 



Here p, f a (p, in) are the incoming momentum and occupancy. 6 Eq. (I3.5P is what we would get 
by performing the z-integration in Eq. (13. 3p . keeping f a constant along a classical trajectory 



5 In general no excitation should be treated as perfectly massless. Therefore, for 7 ^> T/m, the statements 
here have only exponentially small corrections, ~ exp(— jm/T), as opposed to the power (phase space) 
suppressed corrections one would expect if m = strictly. 

6 The integration includes backwards-pointing p z < but they will have exponentially small f(j>) and so 
the error from their inclusion is negligible; the same applies for p z which are positive but so small that 
the particle would reflect. 
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and taking advantage of the absence of reflections. 

Since the wall is at rest and uniform in the transverse directions, both energy and p± are 
conserved. Together with 7 1, this makes evaluating (p z ^ n — p z ,out) simple: 

P 2 z,in + m2 (hi) = pl mt + m 2 (h 2 ), 

P Z M + ^ + 0(mVp 3 ) = p 2 , out + + 0(m V) , 

^P z ,in ^Pz,ont 

m 2 (h 2 ) - m 2 ^) 2 3 2 3 

Pz,m-Pz,out ^ h U{m ft, ,p ± /E ) , (3.oJ 

where in the last line we have used E~ l = —p^ 1 + O (m 2 / E 3 , p\j E 3 ) . The terms ~ E~ 3 will 
be 7 -2 suppressed and can be dropped. Therefore the pressure difference is 

T = E( m '(^) - m lM) I T^S in ) + ^(VT 2 ) • (3-7) 

The integration measure and occupancies are covariant 7 and can be performed in any frame; 
it is most convenient to do so back in the plasma frame. This shows that Eq. (13. 7p is identical 
to Eq. (E2D. 

Therefore the backwards pressure on the interface in the limit 7 ^> 1 is found by replacing 
the thermal effective potential V = V VSiC + V T with V = V vaiC + V^fmean field]. This is the 
main result of this section. 8 

Let us comment quickly on scaling. The density of particles increases as 7, due to Lorentz 
contraction. But the mean momentum that each particle induces, in climbing the wall, goes 
as I/7 (since 5p z ~ m 2 /E ~ I/7), which explains why there is a finite large-7 limit. 

As a check, we can repeat the calculation in the plasma frame. Since the wall passes at 
the speed of light and purely in the +z direction, it leaves each particle's p± and E — p z 
unchanged. 9 However E 2 —p 2 must change by Am 2 , which uniquely determines the changes 
to a particle's 4-momentum: 

. m 2 (/i 2 ) — m 2 (hi) m 2 (ho) — m 2 (hi) \ 

(E,p z ,p ± )^(E + ^-^ ^l 7Pz + ^lL • 3.8 

V 2(E-p z ) 2{E-p z ) J 

In computing the momentum transfer to the phase interface, we must remember to consider 
the flux of particles through the wall, not the particle density. The flux differs from the 



7 Note that p, f(p) are defined based on their symmetric phase [z > 0) values. The expression is covariant 
only if we identify E = E p h i a since that is the energy which satisfies the mass shell condition in the 
symmetric phase, where p is defined. 

8 Some readers might worry that this leading-order analysis may receive large higher-loop effects of form 
a In 7. In particular, since the wall "sees" particles with extremely high energies, shouldn't it analyze 
their partonic content? We believe that the answer is "no." Partonic content becomes important when a 
particle is analyzed not with a large energy, but with a large available transverse momentum needed to 
put the partonic contents on-shell. The wall is uniform in the transverse direction and so cannot impart 
any transverse momentum; so the actual transverse analysis scale is infrared and partonic content is not 
probed. 

9 E — p z conservation is the same as energy conservation in the wall rest frame. 
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particle density by a factor of Relative = 1 — v z = l—p z /E. Therefore the force per unit area 
is 

A - ( 2vr )3M^in^ re lative 2 {E - Pz ) 

which is the same as Eq. (13. 7\i . 



IV. EXAMPLE: SINGLET EXTENSION OF THE STANDARD MODEL 

Here we apply the above results to the simplest extension of the Standard Model which 
can provide a strong electroweak phase transition; the model extended by one real singlet 
scalar field s. This field need not have an s <-> — s discrete symmetry If not, then s 3 and 
sH^H potential terms can appear at tree level, potentially leading to a strongly first order 
electroweak phase transition, as recently studied in 32). 10 

So consider the scalar sector of an extension to the Standard Model; besides the complex 
doublet Higgs field H there is a real singlet field s, with tree-level potential 

V = V H + V HS + V S (4.1) 

with 

V H = -ft 2 h H^H + X h (H^H) 2 (4.2) 

V H8 = \&Hs + ^Hs 2 (4.3) 

Vs = zK s z + ^ + ^ s 4 (4.4) 
2 3 4 K J 

We want to probe whether or not the bubble wall runs away in this model as a function of 

the parameters A, A^, A s , [i 2 ., n 2 h1 a, b, and in particular we want to see how the wall velocity 

is correlated with the strength of the transition (and possibly other attributes). 

To do so we perform a Monte-Carlo sampling over the values of the 7 scalar Lagrangian 

parameters, fixing fi^/Xh to give the right physical Higgs VEV. For each choice of scalar 

parameters, we take a series of steps: 

1. First we must make sure that the parameters give a vacuum consistent with phe- 
nomenology. We enforce absolute stability of the theory, requiring A^ > 0, X s > 0, 
and A > —y/AXfJ^. Then we require that the Higgs vacuum is deeper than any s-only 
vacua and that the scalar states in the Higgs minimum obey experimental Higgs mass 



constraints 34 11 . 



10 

s <-> — s discrete symmetry, the s particles are an interes 
11 We impose a simplified version of the limits found in 



Actually the transition can be strong even without these cubic terms; furthermore, in the theory with 

-inp; dark matter candidate [33| . 
34(. First one computes the masses of the two 
scalar states and their mixing angles with the pure Higgs direction. Then we approximate the constraint 
found in 



34| to require cos 2 9 < 0.25 if to < 115GeV and cos 2 8 < 0.04 if m < 90GeV. 



10 



2. We compute the thermal effective potential and place the following constraints on the 
minimum structure as a function of temperature. The h minimum (by which we mean 
any minimum with a nonzero condensate of the Higgs field) must become unstable 
(become an inflection point rather than a minimum) at some point as the temperature 
is raised; we call this temperature where it suffers spinodal instability T spin i. Then 
we require that the s-minimum which is stable above T sp i n i should become spinodally 
unstable as the temperature is lowered, at some T spin2 > 0. The latter condition 
ensures that the universe does not become trapped in an electroweak-symmetric phase 
as it cools. (The criterion is a little too strong; we really only need ensure that 
the nucleation action to leave the s-phase for the /i-phase is small enough at some 
temperature to allow the transition to occur. Spinodal instability means the nucleation 
rate diverges, so it is a sufficient condition for this to occur.) 

We are mostly interested in the behavior around the transition temperature and in 
cases which are marginal between runaway and finite-velocity bubble walls. There- 
fore we expect moderate scalar VEVs h/T ~ 1, which allows the use of the high 
temperature expansion for the thermal potential, which we take to be 0, 35 



1 - = h 2 T 2 I V + 3g' 2 + I2yj p + 24A fe + 2A \ + ^ f2\ + 3A f 



96 / V 24 

(a + b)sT* h*T f3 + o {g2 + gl2)3/2 | 



12 12tt \A 

Here g,g',yt op are the weak, hypercharge, and top quark Yukawa couplings respec- 
tively; all other terms arise from integrating out the scalars. Note that this potential 
makes a few approximations, dropping the effects of Debye screening and the cubic 
terms for scalar fields. It is not our intention to perform a state of the art analysis, 
and the dropped effects are not too important if the transition is relatively strong 
with relatively light scalar masses near T c , which is the case when the bubble wall is 
marginal between finite velocity and runaway with a transition strong enough to be 
physically interesting (h/T > 1). 

Since very weak transitions are physically uninteresting, we discard all transitions with 
^(7spini)/^spini < 0.2 from further analysis. 

We determine the critical temperature, where the s and h phases have equal free energy. 
Then we find the nucleation temperature T nuc , which we take to be the temperature 
where the energy of a critical bubble is S nuc = E/T = 100. 12 The definition of the 
critical bubble energy is standar d [311] and we use the algorithm for multi-field bubble 
finding presented Appendix 2 of [36]. 



12 It is often stated that bubble nucleation occurs when exp(— S nuc ) — (T/H) 4 so there is one bubble 
nucleation per Hubble volume per Hubble time. This is wrong; one should define T nuc — —dS nuc /dt = 
HTdS nuc /dT, which is the time it takes for the nucleation rate to change appreciably; then one demands 
that exp(— S'nuc) = (T/T) 4 so there are enough bubbles to fill the universe in the time it takes for the 
nucleation rate to change appreciably. This guided our choice of 100 (see also Q). 
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FIG. 3: Example of the effective potential V and pseudo-potential V used to determine bubble wall 
runaway, for two sets of parameters which produce the same value of h/T nuc . The left plot is a case where 
the bubble wall remains of finite velocity; in the right case it runs away. Note the different scales on the 
vertical axis; the scalar excitations in the two minima are much lighter for the left case than for the right. 

4. Using the potential at T nuc and the procedure of the last section, we determine whether 
or not the bubble wall can run away. In detail, we must determine two "potentials" 
at T nuc ; the actual potential V(T nnc ) = V va-C + V T (T nnc ) and the mean-field potential 
expanded about the s minimum using Eq. (14.5j) . V. For our case and thanks to our 
neglect of cubic terms arising from the scalar fields, V(h, s) is the same as Eq. (14.5ft 
without the h 3 T/12ir type terms. 

We plot both the true and mean-field potential for two illustrative cases in Fig. [31 
The black curves are the actual thermal (equilibrium) effective potentials, evaluated 
along the curve in the h, s plane taken by the critical bubble profile. 13 Naturally, in 
both cases they indicate that the /i-phase is deeper (preferred). The red curves are 
the mean-field potentials V defined in the last section, which we can think of as non- 
equilibrium effective potentials. The value of the red curve at a given value of h, s is 
the free energy cost per unit volume to force a phase interface to sweep through the 
plasma, if the phase interface changes the plasma from the h — phase to a phase with 
the given h, s value. For values of (h, s) where it is positive, it would cost energy to 
create a phase with these VEVs if the phase is to be produced by an ultra-relativistic 
interface. Where the potential is negative, there is leftover free energy available from 
creating the /i-phase, which goes into accelerating the bubble wall. Therefore the 
left figure represents a case where the bubble wall cannot run away. The right figure 
represents a case where it will run away. Note however that the value of h where the 
red curve is minimized is not the same as the equilibrium h VEV. This means that, 



V is evaluated along the path through the (h, s) plane which is actually explored by the critical bubble; 
the i-axis is the affine distance in the (h, s) plane along this curve. The cusp in the left figure is at 
the (h, s) value at the center of the critical bubble; we used a straight-line extension from there to the 
minimum of the potential, so the path in the {h, s) plane had a cusp at this point. 
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FIG. 4: Scatter-plots of h/T versus the free energy available in the transition (left) and versus the lighter 
physical vacuum scalar mass (right). In both figures, red dots are parameter values which give runaway 
bubble wall velocities, while black values give finite bubble wall velocities. For the transition to provide a 
large free energy or to arise from relatively heavy scalars, the bubble wall must generally be of the runaway 
type. 



immediately behind the bubble wall, the value of h would be slightly lower than in 
equilibrium. The VEV h will relax to the equilibrium value as the non-equilibrium 
state after the bubble wall passage re-equilibrates. The height indicated as A in the 
diagram is free energy density available to accelerating the phase interface; the height 
indicated as B is free energy density lost to entropy production in re-equilibrating the 
system after the wall's passage. 



Some results from our Monte-Carlo study of the parameter space of this toy model are 
shown in Figure HI The study shows, as expected, that weak transitions generally involve 
finite wall velocities, while strong transitions almost always require runaway bubble wall 
velocities. However a small subset of strong transitions (as measured by h/T, the relevant 
quantity for baryogenesis since it determines whether baryon number is violated after the 
transition 37|) have finite bubble wall velocities rather than runaways. These strong but 
non-runaway transitions occur when the potential between the two (s and h) minima is very 
flat. The figure also shows that, in these cases, there is also generally a light physical scalar 
state in vacuum, which must be nearly aligned with the s direction in order to meet Higgs 
search bounds. (The cutoff in h/T at small values is imposed by our procedure since we 
discard very weak transitions.) 



V. CONCLUSIONS 



The physical consequences of the electroweak phase transition depend on the propagation 
speed of the phase interface. Theories with singlet scalars can have very strong transitions 
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without requiring large couplings or violating the Higgs mass bound, because the potential 
can already possess multiple metastable minima at the mean-field level. We show here that 
such strong mean-field type transitions generally - but not always - involve interfaces which 
propagate at virtually the speed of light, as occurs in vacuum transitions. However, in the 
MSSM or other theories where the first order nature of the transition is fluctuation induced, 
the bubble wall can never run away. We provide a clean, simple criterion, valid at the 1-loop 
level of analysis, which determines whether or not bubble walls run away. Specifically, one 
defines V, the vacuum effective potential V V3JC plus the quadratic-order Taylor expansion 
of the thermal part V T of the potential about the high-temperature phase. If the low- 
temperature phase is preferred at T nuc by V, then the wall runs away; otherwise it has a 
finite terminal velocity. For the case of a real singlet extension of the Standard Model, we 
find that very strong phase transitions (say, h/T > 2) either have runaway bubble walls or 
have very slight supercooling, a small available free energy at the transition, and very light 
physical scalar excitations. 

It would be interesting to apply this work to the nMSSM. For cases where couplings are 
relatively large or where two- loop thermal effects are important, it would also be necessary 
to understand how to incorporate these effects in to our treatment. 
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APPENDIX A: THERMODYNAMICS OF RUNAWAY WALLS 

This Appendix explores the thermodynamics which establishes whether runaway bubble 
walls are thermodynamically consistent. Consider two phases, with pressures P\{T) and 
P^iry, the entropy density s = dP/dT and energy density e = sT — P follow from P(T). 
Here P\{T) is for the high temperature s-phase, p2(T) for the low-temperature h- phase, so 
for T < T c we have P2(T) > P\{T). The phases are separated by a planar interface moving 
in the +z direction at the speed of light, with surface energy density a(t). We work in the 
rest frame of the unburnt s-phase plasma; the /i-phase plasma fills the region z < t and flows 
with forward velocity v 2 - The temperature of the unburnt plasma is Ti, the burnt plasma is 
T 2 . We want to determine for what values of Ti the interface gains surface energy density, 
da/dt > 0. 

This problem has three unknowns; T 2 , v 2 , and da/dt. We have two equations and one 
inequality with which to determine them; conservation of energy and of z-momentum, and 
the positivity of entropy production. The equation of state provides the stress tensor in each 
phase: 



e, = €i(Ti), Pi = PiiTi) . 



u x = (1, 0, 0, 0), u 2 = (72, 0, 0, 72^2), 



(Al) 
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The index i indicates which phase's equation of state is to be used. Per unit time, the energy 
per area of phase 2 increases by T 2 °, while the energy per area of phase 1 decreases by T®°. 
These must be balanced by the momentum flows and the energy change of the bubble wall: 

— - l x 1 2 +1 2 1, 

= ei - (72 - ljv 2 )e 2 + (1 - 7 2 + ^v 2 )P 2 ■ (A2) 
Similarly the z-momentum must be balanced: 

rpzO rpzO I rpzz rpzz 

= -P 1 + (7 2 \ 2 - 7 2 t; 2 )e 2 + (1 + - ^v 2 )P 2 . (A3) 
Adding these yields a simpler expression: 

1% = e 1 - P x - e 2 + P 2 . (A4) 
at 

The total entropy of the system must also increase. The entropy density of phase 2 is 
72(^2 + P-i)jTi (the factor 72 because of Lorentz contraction) and the volume increase of 
phase 2 per unit time is (1 — v 2 ). Therefore entropy increase requires 

l2{1 - V2 f^t P ^ > ^t P ^ (A5) 

where the RHS is the reduction of entropy stored by the shrinking s-phase. 

The bubble wall runaway is permitted if there is a solution with da/dt > 0, so energy 
accumulates in the bubble wall rather than being lost. Unfortunately two equations and an 
inequality are not enough to solve for the three unknowns da/dt, v 2 , and T 2 . If we somehow 
knew how much entropy is generated at the wall, we could replace Eq. (IA5I) with an equation 
and all variables would be determined. This requires some microscopic dynamics details, 
which are treated in the main text of the paper. Here we will just determine under what 
conditions a runaway wall is thermodynamically consistent. That is, we can determine the 
domain of T\ values where it is at least thermodynamically consistent for a bubble wall 
to run away. We do this by setting the entropy production to zero, which determines the 
bubble wall velocity: treated as an equality, Eq. (IA5I) becomes 



1-x 2 _ (e 1 + P 1 )T 2 Sl 

1+x 2 (e 2 + P 2 )Ti s 2 

Then the remaining two equations determine da/dt in terms of the initial temperature I\. 

Specifically, since the /i-phase becomes more favorable as the temperature drops, we 
can search for the temperature at which ultra-relativistic bubble propagation first becomes 
thermodynamically feasible, which is where da/dt = at zero entropy generation. Then T 2 
is determined from T\ by using Eq. (1A4I) : 62(T 2 ) — P 2 (T 2 ) = ei(Ti) — Pi(Ti). The critical 
value of T\ is the value where the so-determined T 2 , v 2 will give a valid solution to Eq. (lA2h . 

We evaluated this value T% for each set of parameter values in the toy model discussed in 
Sec. [TV] We also determined T reheat , determined by the criterion 

ei(T rehea t) = e 2 (T c ) . (A7) 
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That is, T re heat is the temperature where the supercooled s-phase has the same energy density 
as the /i-phase at T c . Its physical importance is that, for slow bubble walls, if T nuc > T re heat 
then the universe will heat back up to T c before the phase transition completes and there will 
be a period of mixed phase; whereas if T nuc < T reheat then the supercooling is enough to use 
up all of the latent heat of the transition without reheating to T c and no quasi-equilibrium 
mixed-phase situation occurs. 

Surprisingly, we find in the toy model that in every case T c > T\ > T re h ea t- That is, the 
criterion that ultra-relativistic bubble walls are thermodynamically consistent is in practice 
a weaker condition than the condition of sufficient supercooling not to reheat to T c . Further, 
we find T reheat > T nuc for all but a very small corner of the parameter space of the theory. 
Therefore, in the singlet-scalar model, supercooling proceeded below either temperature 
throughout almost the entire parameter space which provides a first order phase transition. 
The fact that many parameters give bubble walls which cannot run away means that the 
production of entropy at the bubble wall significantly reduces the range in which walls run 
away, relative to what is thermodynamically consistent. 
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